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^- Abstract. Consider random shadows of a cube and of a regular tetrahe- 

dron. Area and perimeter of the former are positively dependent (with correla- 
tion 0.915...), whereas area and perimeter of the latter appear to be negatively 
psj ■ dependent. This is only one result of many, all involving generalizations of 

mean width. 

^ ■ Let C be a convex body in M". Let S" be a l-dimensional subspace passing through 

the origin. A width is the length of the orthogonal projection of C on S. If S is 
uniformly distributed on the Grassmannian manifold G"'^, then the width w{S) is a 
random variable. First and second moments of w are known for C = the regular 
ra-simplex, C = the ra-cube and C = the regular ra-crosspolytope, each centered at the 
origin [Il|2]. The mean width is sometimes called the mean linear projection [3]- 
a slight abuse of language - because a projection of C on 5* has only one parameter 
^ ! of interest: its length. 

^ I Let instead S* be a 2-dimensional subspace passing through the origin. Orthogonal 

Tjj- ' projections of C on S" now have three parameters of interest. One parameter is the 

CN ■ number of polygonal vertices; associated probabilities will be mentioned later. The 

cn . other two parameters constitute direct generalizations of width and will be our main 

^ ! focus: 

• a chorowidth is the polygonal area of the orthogonal projection of C on 5 
k>( ■ ( "choro" is Greek for place or area) ; 

c^ ■ • a periwidth is the polygonal circumference of the orthogonal projection of C 

on 5* ( "peri" is Greek for around or enclosing) . 

If S is uniformly distributed on G"'^, then chorowidth cw{S) and periwidth pw{S) are 
random variables. We shall examine joint moments of {cw,pw) for the same regular 
polytopes as before, for the special cases n = 3 and n = 4. 
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1. Method for Computing Moments 
To generate a random 2-subspace S* in R^, we select a random point U uniformly on 
the 2-sphere of unit radius. The desired subspace is the set of all vectors orthogonal 
toU. 

We then project the (fixed) convex polyhedron C orthogonally onto S. This is 
done by forming the convex hull of images of all vertices of C. If C = the regular 
3-simplex, the resultant polygon in the plane has 3 or 4 vertices almost surely. If C = 
the 3-cube, there are 6 vertices almost surely; if C = the regular 3-crosspolytope, 
there are 4 or 6 vertices almost surely. We compute the area and circumference of 
this polygon via standard formulas. 

More precisely, ii U = {x, y, z) is of unit length, then the matrix 



M, 




xy 



y 



yjl—x'^ \/l — x'^ 





projects C orthogonally onto a plane, rotated in R^ to coincide with the 2-subspace 
spanned by (1, 0, 0) and (0, 1, 0) for convenience. Let T be the 3-row matrix whose 
columns constitute all vertices of C . Then the first 2 rows of M^T constitute all images 
of the vertices in R^ and 2-dimensional convex hull algorithms apply naturally. 
Spherical coordinates in R^: 

X = cos 9 sin if , y = sin^sin(/9, z = cosip 

will be used throughout for f/, where < ^ < 27r, < (/? < vr. The corresponding 
Jacobian determinant is sin ip\ it is best to think of (^, if) as possessing joint density 

Moving up a dimension, to generate a random 2-subspace S* in R^, we first select 
a random point U uniformly on the 3-sphere of unit radius. The set of unit vectors 
orthogonal to U form a 2-sphere; we select a random point V uniformly on this 2- 
sphere. The desired subspace is then the set of all vectors orthogonal to both U and 
V. 

We next project the (fixed) convex polyhedron C orthogonally onto S. This is 
done by forming the convex hull of images of all vertices of C. If C = the regular 
4-simplex, the resultant polygon in the plane has 3, 4 or 5 vertices almost surely. If 
C = the 4-cube, there are 8 vertices almost surely; if C = the regular 4-crosspolytope, 
there are 4, 6 or 8 vertices almost surely. Area and circumference of this polygon are 
computed as before. 

More precisely, ii U = {x, y, z, w) and V = (p, q, r, s) are orthogonal and of unit 
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length, then the matrix 



Ma 



/^/I 



V 



p^ — x"^ 








pq+xy 



pr+x z 



rw—sz 
■UX—p^—x^ 






a/1— p2— x^ 
qw—sy 

Wl—p'^—x'^ 






ps+xw \ 
q z~ry 



■i/l—p'^—x^ 



Wl—p'^—x'^ 






/ 



projects C orthogonally onto a plane, rotated in R^ to coincide with the 2-subspace 
spanned by (1,0,0,0) and (0,1,0,0) for convenience. Let T be the 4-row matrix 
whose columns constitute all vertices of C. Then the first 2 rows of M^T constitute 
all images of the vertices in M^ and 2-dimensional convex hull algorithms again apply 
naturally. 

Spherical coordinates in W^: 

X = cos 6 simp sin ip , y = sin 6 sirup sin ip, z = cos ip sin ip , w = cosip 

will be used throughout for U, where < ^ < 27r, 0<(p<7T,0<ilj<7T. It follows 
that V is given by 



(p\ 




( "M 




/ -z\ 




( -w\ 


r 


= COS K, sin A 


X 


+ sin K, sin A 


w 

X 


+ cos A 


y 


V V 




\ ^ / 




\-y / 




\ X J 



where < /t < 27r, < A < vr. It is best to think of {9, ip, ip, n, A) as possessing joint 
density ^ sin ip sin^ ip^sinX. 

2. Three-Dimension AL Results 
From the foregoing, we calculate mean chorowidth and mean periwidth via integra- 
tion: 

E (cw) = area (convex_hull (M3T)) — sin (p dip dO, 

J J 47r 







2-7r TT 



E [pw) = circumference (convex_hull (M3T)) — sin (p dip dO 


and likewise for mean square chorowidth, mean square periwidth and joint moment. 
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2.1. 3-Simplex (Tetrahedron). 

E{cw) = ^= 0.433012701892219... = !!^!l^^^^!^, 

^ ^ 4 4 

1 /9 1 

E [cw^) = -^ + j ^ arcsec(3) = 0.188561220515812..., 

E (pw) = I (vr - arcsec(3)) = 2.865949854373527..., 

E{pw^) =8.2170808733..., 

E {cw ■ pw) = 1.2406348222... 

which imply that the correlation between cw and pw is ~ —0.188. Closed-form 
expressions for E {pw^) and E {cw ■ pw) are not known. The fact that 

surface area of a \ / mean area of 

3D convex body J y its 2D shadow 

was first noticed by Cauchy [U El El [7] . We also have [H [9] 

2 
P (projection has 3 vertices) = — (3 arcsec (3) — tt) = 0.3509593121..., 

n 

3 

P (projection has 4 vertices) = - (tt - 2 arcsec (3)) = 0.6490406878.... 

n 

2.2. 3-Cube. 

^ , , 3 surface area 

E (cw) = - = 1.5 = , 

^ ' 2 4 ' 

E (cw^) = 1 + - = 2.273239544735162..., 

n 

E{pw) = -TT = 4.712388980384689..., 

E (pw^) = 8 + 6713^2 (-|, i, |; 1, 2; l) = 22.237117433439470..., 
E{cw-pw) = 2 + — = 7.092958178940650... 

TC 
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which imply that the correlation between cw and pw is 0.915.... We shall present 
more details underlying these results later, especially that involving the generalized 
hjrpergeometric function 

,F,ia„ a„ as; h, b,; z) - r(^^)r(a,)r(a3) ^ T{h, + k)T{h, + k) B/ 

whose appearance is quite unexpected. 
2.3. 3-Crosspolytope (Octahedron). 

E(c^) = ^ = 0.866025403784438... = '^'^^'^^ ^'^^ 

1 /9 1 

E icw'^) = - + - — arcsec(3) = 0.754244882063249..., 

2 7r 27r 

¥.{pw) = 3arcsec(3) = 3.692878252022324..., 

E(pw;2) = 13.6639421274..., 

E {cw ■ pw) = 3.2074623048... 

which imply that the correlation between cw and pw is ~ 0.878. The fact that 
cws-crosspoiytope behaves like 2 • cw3_simpiex will be discussed shortly. We also have [10] 

3 

P (projection has 4 vertices) = - (vr - 2 arcsec (3)) = 0.6490406878..., 

TT 

2 

P (projection has 6 vertices) = — (3 arcsec (3) — vr) = 0.3509593121.... 

vr 

See [m US |T3] for related probabihties governing planar cross sections of the tetra- 
hedron, cube and octahedron (rather than projections). 

3. Tables of Intrinsic Volumes 
Let n be a rectangular 4-parallelepiped in R^ of dimensions Zi, Z2, z^, Z/^. It is 
well-known that [H] 

^4(0) = ZxZ2Z'iZ^, 

V3(n) = ZxZ2Z-i + ZxZ2Z^ + Z^Z-^Z^ + 2:2^3^4, 

V2(n) = ZxZ2 + ZxZ-i + 2:1^4 + 22-^3 + -22-24 + -23^4, 

Vl(n) =2:1+2:2 + 2:3 + ^4 
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are the elementary symmetric polynomials in four variables. In M", there are n 
such intrinsic volumes, corresponding to the n elementary symmetric polynomials. 
Limiting approximation arguments enable us to compute Vj{C) for arbitrary convex 
C. 

One motivation for our work is to generalize the following two tables [1] for n = 2: 



and n = 3: 



to n = 4: 



V2 = area, 2Vi = circumference 



V3 = volume, 2V2 = surface area, ^Vi = mean width 



V4 = hypervolume, 

2V3 = hyper- surface area, 

|V2 = mean chorowidth, 

|Vi = mean periwidth = vr (mean width) . 

Another motivation is to use formulas for V2 in [151 Il3| to deduce that 



00 „ 

n--2 



^, , n{n + l) f _3^2 /l + erf(x) . , 

E (cw„_simpiex) = \, ^ r \ 2^ ' 



—00 

for the regular n-simplex in R"; 

n 
E (cw;„_cube) = - 

for the n-cube in M"; and 

lUi l^Cti'jj-crosspolytope j ^ / 6 Cri yX J dX 



for the regular n-crosspolytope in M". In contrast, E (pWn) are obtained simply by 
forming the product vr ■ E(w„), and associated E (w„) are tabulated in [H |2]. No 
such general expressions are available for higher moments. 
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4. Four-Dimensional Results 
From the foregoing, we calculate mean chorowidth and mean periwidth via integra- 
tion: 

27r 7r 27r n n 

K{cw) = / / / / / area (convex_hull (M4T)) — - sin if siir^ ^p sin X d^p d^p d9 dX dK, 
00000 

27r TV 2-77 7r tt 

E (pw) = / / / / / circumference (convex_hull (M4T)) — ^ sin </) sin^ V^ sin A c??/^ ^99 d^ o?A rfK 

00000 
and likewise for mean square chorowidth, mean square periwidth and joint moment. 

4.1. 4-Simplex. 

5^/3 

E{cw) = -i— (tt - arcsec(4)) = 0.418889720727840..., 

E{cw^) =0.176..., 

E{pw) = — (27r - 3arcsec(3)) = 2.748401146360593..., 
37r 

E{pw^) =7.56..., 
E {cw ■ pw) = 1.15... 

which imply that the correlation between cw and pw is ~ 0.1. Closed-form expressions 

for E{cw^), E [pw"^) and K{cw ■ pw) are not known. We also have 

P (projection has 3 vertices) ~ 0.146, 

P (projection has 4 vertices) ~ 0.585, 

P (projection has 5 vertices) ~ 0.269. 

These values are consistent with a theorem in |H] that the expected number of vertices 
should be 

I — °° 3 

'°v!/^""^ C^f^"0 ^" = 10(l-Aarcsec(3))=4.122... 
—00 

^ 3(0.146) + 4(0.585) + 5(0.269). 
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4.2. 4-Cube. 



E (cw) = 2 



E{cw^) =4.04..., 

E{pw) = - = 5.3, 

E{pw^) =28.4..., 
E {cw ■ pw) = 10.7... 
which imply that the correlation between cw and pw is ~ 0.9. 

4.3. 4-Crosspolytope. 

E(cw) = -^ = 0.769800358919501..., 

E{cw^) =0.598..., 

16 

E{pw) = —{71 - 2arcsec(3)) = 3.461550335020567..., 

n 

E{pw'^) = 12.0..., 
E {cw ■ pw) = 2.67... 
which imply that the correlation between cw and pw is ~ 0.8. We also have 

P (projection has 4 vertices) ~ 0.463, 

P (projection has 6 vertices) ~ 0.478, 

P (projection has 8 vertices) ~ 0.059. 

These values are consistent with a theorem in [10] that the expected number of vertices 
should be 

oo 

2 f _22;2 „, .0 , ^ . / . 2 



48^/- / e-^"" eTi{xydx = 24 ( 1 - -arcsec(3) ) = 5.192... 





4(0.463) + 6(0.478) + 8(0.059). 
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5. Further Work 
After having written the preceding, we discovered [17], which provides new insights in 
the 3-diniensional case. Earher papers in this hne of thought include [HI [191 EOl [21] . 
Consider the 3-cube with vertices 

(1 _i 1 1 _i _i 1 _i 
2222 222 2 

1 1 _1 1 _1 1 _1 _1 

22 22 22 22 

111 _1 1 _1 _1 _1 

222 22 222 

A starting point here is the simuhaneous system of equations 

cw = X + y + z 

PW = 2 (Vl - X^ + y/1 - y2 + VI - zA 

given a unit vector U = {x, y, z) in the first octant. This makes possible, for example, 
the derivation of a closed-form marginal density for cw (although not a joint density). 
Let us focus on computing E (pw^) and E {cw ■ pw). Contributing to E (pw^) are three 
terms like 

tt/2 tt/2 

IT 



/ (l — cos^ 9 sin^ (/?) sin ip dO d^p 



I . . ,_ .,--. -^- -^. 

3 





and six terms like 

7r/2 7r/2 

J = y 1 — siii^ 9 sm"^ (f y 1 — cos"^ if sin(f d9 dip 



tt/2 



/2 
E(sin(^)sinVciv' = y3i^2(-|,i|;l,2;l) 





where 



7r/2 1 , 

r2+2 



E{0 = J^l-eM0rd9 = j^\^dt 



is the complete elliptic integral of the second kind. The final result is 32(3/+6 J)/(47r) 
For E {cw ■ pw), the calculations are simpler, with / = 7r/6 and J = 2/3. 
An analogous system of equations for the regular 3-simplex with vertices 

^3 6 6 

r= I I -i 

4 12 12 12 
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is more difficult. Define constants 

7 = arccot (2^2) = 0.339..., 6 = arccot (^2) = 0.615... 

and functions 

, X f arcsec ('2A/2tan(/9) if 7 < (/? < 7r/2, „, . 27r 

"(^) = lo ifo<v.<7, /3(v^) = Y-«M- 

Given a unit vector U = {x, y, z) in tlie ffist dodecant (one of twelve regions in 
3-space), we have 



12 
cw 



j^ {VQx + 3^21/ + 73^) if (5 < V? < f and (5{^) < ^ < f , 

^ (y2a; + z) if 7 < V? < f and < ^ < min{a(v9), /3(v?)}, 



^z if < ^ < (5 and ^(v?) < ^ < f 

and the three respective expressions for pw are 

\ I 3W4 - (V^x - yY + V3J12 - (x- VSy + 2^2^^ + 273^3 - (x - v^zV ] 
^ I VsJa - (^VSx -yy + VsJa - (ySx + yy+ 



12 



(x - VSy + 2\f2z\ + W 12 - [x \ VSy + 2V2z) j 



^ ^ -U- (V3x-y] +2x/l^^+\U- (V3x + 



An analogous system of equations for the regular 3-crosspolytope with vertices 



T 



/^ -^ \ 

4s — ^ 



V u u u u 73 -73 y 

is similar. Define functions 

, X f 7r/4 — arcsec fv/2tan(/9) if 5 < (/? < 7r/2, „. x . x 

"(^) = lvr/4 if0<v.<5, /^(V^) = -"(V^) 



2 . 
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where the constant S is the same as before. Given a unit vector U = {x, y, z) in the 
first hexadecant (one of sixteen regions in 3-space), we have 

( \{x + y + z) if (5 < (/? < f and max{a((/?), /^(y?)} < ^ < f , 
cw= I X iff<V^<fandO<^< /3((^), 

\ z ifO<(^<fandO<e< a(^) 

and the three respective expressions for pw are 

V2 ( V2 - (x + i/)2 + v/2 - (x + ^)2 + v/2 - (y + z) 

^i^^2-{y-zY+^2-{y + zf), 

V2(^^2- {x -yf + ^2- {x + yf) . 

A rigorous proof that cw3_crosspoiytope and 2 ■ cu73_simpiex are identically distributed re- 
mains open. 

Let us return finally to (ordinary) width w and questions left unanswered in 
The 2-cube (square) and regular 2-simplex (equilateral triangle) have vertices 



3 

respectively. Given a unit vector (x, y) in the first quadrant, w = x + y for the square. 
Given a unit vector {x,y) in the first sextant, w = ^x + ^y for the triangle. The 
corresponding width densities are 

-^=i= ifl<^<y2, j ^^^= ii^<w<l, 

otherwise, y otherwise 

respectively, via simple argument. 
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